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Ijr^ ' Abstract. The discovery of the X-ray source IGR J17252-3616 by INTEGRAL was reported on 9 February 2004. Regular 

monitoring by INTEGRAL shows that IGR J17252-3616 is a persistent hard X-ray source with an average count rate of 0.96 
counts s"' (~6.4 mCrab) in the 20-60 keV energy band. A follow-up observation with XMM-Newton which was performed on 
21 March 21 2004, showed that the source is located at R.A. (2000.0) = 1 7'' 25"' 1 1.4' and Dec. = -36°16'58.6" with an uncer- 
tainty of 4". The only infra-red counterpart to be found within Ihe XMM-Newton error circle was 2MASS J17251 139-3616575, 
which has a Ks-band magnitude of 10.7 and is located 1" away from the XMM-Newton position. 

The analysis of the combined INTEGRAL and XMM-Newton observations shows that the source is a binary X-ray pulsar with 
a spin period of 413.7 s and an orbital period of 9.72 days. The spectrum can be fitted with a flat power law plus an energy cut 
off (F ~ 0.02, £c ~ 8.2 keV) or a Comptonized model {kT^ ~ 5.5 keV, t ~ 7.8). The spectrum also indicates a large hydrogen 
column density of A'h ~ 15 10~ atoms cm"^ suggesting an intrinsic absorption. The Fe Ka line at 6.4 keV is clearly detected. 
< I ■ Phase-resolved spectroscopy does not show any variation in the continuum except the total emitted flux. The absorption is 

' constant along the pulse phase. This source can be associated with EXO 1722-363 as both systems show common timing and 

, spectral features. The observations suggest that the source is a wind-fed accreting pulsar accompanied by a supergiant star. 
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1 . Introduction ject. The orbits of Be binary systems are eccentric and gener- 
ally have a longer period than OB supergiant systems 15 

X-ray binaries consist of a compact object, either a black hole ^^^^y r^^^^ x-ray emission of HMXB presents a wide variety 

or a neutron star, accreting matter from a companion star; they patterns: from transient to persistent, outbursts on differ- 

are usually classified according to the mass of the compan- ^-^^^ ^^^j^^ (seconds, days and/or years), periodic modu- 

lon as a high mass (HMXB, Mc ^ 10 Mq), intermediate- lations, eclipses and others. Strongly magnetized neutron stars 

mass (IMXB, Mc - 1-10 M^) or low mass X-ray binary accompanied by massive stars show periodic pulsations; the ac- 

(LMXB, Mc < 1 Mo). The LMXB and IMXB accrete mat- ^^^^^^ ^^^^^^ funneled towards the poles by the magnetic 

ter through Roche-Lobe overflow from the companion star ^^j^ j^^^j^g ^„ ^^^^^^^^ observed X-ray emission 

and through an accretion disk around the compact object ^j^^^^ ^^^-^^^ ^^^^^ ^^e line of sight. These HMXB usu- 

(Tauns & van den Heuvel 2005). ^jjy ^j^^^ ^ j^^j.^ ^-ray spectrum between 2-10 keV with an 

The HMXB can be divided into two categories as the com- energy cutoff around 10 keV. 
panion star can be either an OB supergiant or a Be star. OB 

supergiant stars feed the compact object through strong stel- Several new hard X-ray sources have been discovered 

lar winds and/or in some cases Roche-Lobe overflow. Be stars by INTEGRAL (Winkler el al. 2003) in surveys of the galac- 

expell matter around their equator that fuels the compact ob- tic plane. Most of them have been detected by IBIS/ISGRI 

(Ubertini et al. 2003; Lebrun et al. 2003), the most sensitive 

Send offprint requests to: I.A. Zurita Heras instrument on board INTEGRAL between 20 and 300 keV. 

e-mail: Juan. ZuritaSobs.unige.ch The brightest sources detected during the first year of the 
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INTEGRAL mission are listed in Bird et al. (2004). A few tens 
of them have never been detected before INTEGRAL'S ob- 
servation. Most of these new objects show common features 
in their spectra, such as a high intrinsic low-energy absorp- 
tion, and they are believed to be HMXB (Walter et al. 2003; 
Rodriguez et al. 2003; Patel et al. 2004). 

IGR J 1 7252-36 16's discovery was reported on February 
9, 2004, with 13 other hard X-ray sources (Walter et al. 
2004) detected with ISGRI in all-sky mosaic images built 
from core programme data. However, there is evidence that 
IGR J17252-3616 has already been observed with previous 
missions. A galactic plane scan performed by EXOSAT in June 
1984 revealed a point like X-ray source emission, GPS 1722- 
363 (Warwick et al. 1988), at a position not compatible with 
IGR J 17252-36 16 but with a low accuracy. Ginga observations 
in 1987 and 1988 confirmed the presence of a powerful X-ray 
accretor neutron star, XI 722-36, with a pulsation of 413.9 s and 
important variations of the intensity in X-rays (Tawaraet al. 
1989). The source also showed a hard spectrum with important 
low-energy absorption and an emission line at 6.2 keV. Further 
investigations with Ginga in 1988 confirmed the spectral anal- 
ysis, and lower limits for the orbital period of 9 days and the 
mass of the primary star of 15 Mgwere deduced from pulse 
timing analysis (Takeuchi et al. 1990). Both investigations con- 
clude that the system is a HMXB. Recently, those results were 
confirmed by Corbet et al. (2005) using RXTE data. They also 
found an orbital period of 9.741+0.004 d. 

IGR J17252-3616 is regularly monitored by INTEGRAL. 
A follow-up observation lasting three hours was performed 
with the X-ray Multi-Mirror Mission (XMM-Newton) on 21 
March 2004. The new data available on this source is presented 
in this paper. INTEGRAL and XMM-Newton observations and 
data analysis techniques are described in Sects. 2 and 3, re- 
spectively. The results are presented in Sect. 4 and discussed in 
Sect. 5. The conclusions are presented in Sect. 6. 

2. Observations 

2.1. INTEGRAL 

INTEGRAL is a hard X-ray and y-ray observatory of the 
European Space Agency (ESA) launched on 17 October 2002. 
The payload consists of four instruments: the imager IBIS with 
two detector layers, ISGRI (20 keV-1 MeV) and PICsIT (200 
keV-10 MeV, Labanti et al. (2003)); the spectrometer SPI (20 
keV-8 MeV, Vedrenne et al. (2003)); the X-ray monitor JEM- 
X (3-30 keV, Lund et al. (2003)) and the optical camera OMC 
(V filter, Mas-Hesse et al. (2003)). Most of the INTEGRAL ob- 
serving time is spent in the Galactic Plane. 

IGR J17252-3616 is located close to the Galactic Centre, 
which has been observed regularly. The observing strategy con- 
sists of pointings each lasting ~ 30 minutes distributed in var- 
ious grids around the Galactic Plane. The focus of the effort 
for this paper is on IBIS/ISGRI, since IBIS/PICsIT and SPI are 
less sensitive at energies lower than 300 keV where the source 
is detected. We did not use data from the JEM-X instrument 
because its smaller field of view and the dithering observa- 
tion approach of INTEGRAL mean that the effective exposure 



on IGR J 17252-36 16 is ten times smaller than obtained with 
ISGRI. 

The data set consists of core programme data obtained until 
MJD 53341.1 and of pubUc data obtained until MJD 52928.3, 
giving a total exposure of 6.5 Ms. Table 1 lists the source vis- 
ibility periods, and for each of them the fraction of the time 
when the source was effectively observed in the partially-coded 
field of view (PCFOV). 

2.2. XMM-Newton 

IGR J17252-3616 was observed by XMM-Newton 
(Jansen et al. 2001) on 21 March 2004, from 13:02:45 to 
16:04:45 UTC (MJD 53085.544-53085.671) for a total expo- 
sure of 11 ks. XMM-Newton operates with three instruments 
that cover the optical/UV and X-ray spectral bands. The main 
instrument for the purpose of this work is the European Photon 
Imaging Camera (EPIC) instrument, which consists of two 
MOS (Turner el al. 2001) and one pn (Struder et al. 2001) 
CCD cameras. EPIC has a 30' field of view, and coverage 
of the 0.15-12 keV energy range with imaging, timing and 
spectral capabilities. The EPIC/MOS[12] and pn were all 
operating in imaging science mode with a large window and a 
medium filter 

The XMM-Newton observation was simultaneous with 
three INTEGRAL pointings in revolution 175 between 
MJD 53085.542 and 53085.667. However, the signal to noise 
(S/N) of ISGRI data during those simultaneous observations 
was too low to perform cross analysis. 

3. Data Analysis 

3.1. INTEGRAL 

The data were reduced with the Offline Scientific Analysis 
version 4.2 software (OSA 4.2) that is pubUcly released by 
the INTEGRAL Science Data Centre (ISDC) (Courvoisier et al. 
2003). All the pointings for which the source was either within 
the ISGRI fully coded (FCFOV, 9° x 9°) or partially coded 
(PCFOV, 29° X 29°) fields of view were analysed. Around 
3000 pointings, distributed between revolutions 37 and 244, 
were selected. Sky images were extracted for each point- 
ing and combined into mosaic images with longer exposures. 
As the source is weak, light curves were extracted from the 
imaging results (pointings and mosaics). During flares, light 
curves were extracted on shorter time scales using ii .light 
v7.3. Extracting high-energy spectra with OSA 4.2 remains a 
difficult task for faint sources. Spectra were extracted using 
ii_spectra_extract v2.3.1 and from the mosaic images, and 
compared. In the energy range in which the source is detected, 
both methods give consistent results. The redistribution ma- 
trix and ancillary response files (RMF and ARE) used were 
isgr_rin£_grp_0(S12 . fits and isgr_ar£_rsp_0IS06.fits, 
respectively. The RMF was rebinned in to 17 channels. Average 
fluxes and the source position were obtained from the mosaic 
images with mosaic_spec vl.O'. 

' ii.light and ii_spectra_extract are OSA 4.2 standard tools. 
mosaic_spec is released with OSA 5.0. 
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Table 1. The INTEGRAL Data set. We selected all public and core programme revolutions when IGR J17252-3616 was within 
the ISGRJ PCFOV of 29° x 29°. The observing time fraction was calculated as the ratio between the net exposure time on the 
source and the elapsed time during the visibility periods. 



Visibility Periods 


Start 


End 


Revolutions 


Observing time fraction 




[MJD] 


[MJD] 






1 


52671 


52752 


37- 63 


0.23 


2 


52859 


52925 


100-121 


0.51 


3 


53050 


53116 


164-185 


0.20 


4 


53236 


53294 


226-244 


0.16 



3.2. XMM-Newton 

The Science Analysis System (SAS) version 6.1.0 was used 
to produce new event lists for the EPIC instrument running 
epchain for pn and emchain for MOS[12]. The event Hsts 
were corrected for enhanced background features at energies 
higher than 10 keV, disregarding time lapses when count rates 
above 10 keV exceeded 4 counts s"' for pn, 1 counts"' for 
MOSl and 1 counts ' for MOS2. Finally, only 5.9 ks of 9.2 
ks, 6.2 ks of 10.7 ks and 7.0 ks of 10.7 ks for each instrument, 
respectively, were kept as good time intervals. 

Images were built from the cleaned event lists for the 
MOS[12] and pn cameras with 2" and 4" resolution, respec- 
tively. Bad pixels were disregarded and good events were se- 
lected until the quadruple level. The lower threshold was fixed 
at 0.8 keV as recommended in the calibration status documen- 
tation (see calibration document XMM-SOC-CAL-TN-00 18^ , 
p25). 

To find the source location, the SAS task edetect_chain 
was used on each individual EPIC camera. Four images, with 
energy ranges of 0.5-2, 2^.5, 4.5-7.5 and 7.5-12 keV, were 
created. The source position was calculated as the mean of the 
best position from each MOS and pn camera. 

The source was clearly visible in the pn images in CCD 1, 
near the read-out node. An event list of source+background 
counts was selected from a circle of 50" around the bright 
object. The background was estimated from a region adjacent 
to the source region, with a similar size and at the same dis- 
tance from the read-out node in the same CCD. To build light 
curves, only single and double events with an energy between 
0.4 and 10 keV were collected, as advocated in the user guide. 
The background light curve was subtracted using the FTOOLS 
Icmath. 

Source+background and background spectra were ex- 
tracted from all EPIC cameras, MOS[12] and pn, disregarding 
bad pixels and selecting single and double events. Specific re- 
distribution matrix and ancillary response files were generated 
for each EPIC instrument with the standard tasks rmfgen and 
arf gen, respectively. An average spectrum of the source from 
the complete observation was first obtained and the Xspec ver- 
sion 11.3.1 package was used to fit and plot the resulting spec- 
tra corrected for the background. 

For the phase-resolved spectra, the Xselect version 2.2 
software was used to select events corresponding to a specific 

^ http://xmin.vilspa.esa.es/external/xmm_sw_cal/ 
calib/documentation/index . shtml#EPIC 




Fig. 1. Mosaic image of IGR J17252-3616 with ISGRI during 
revolution 106 (MJD 52877.4-52880.4). 

phase bin given the start epoch and the periodicity of the ob- 
ject. Source and background spectra for each selected phase bin 
were calculated from the same regions as defined previously. 

4. Results 

4.1. Source position 

The source's brightest epoch occurred during revolution 106 
between MJD 52877.4 and 52880.4 with an average flux of 
4.08+0.09 counts s ' and a significance of 47cr. The ISGRI 
20-60 keV mosaic of that epoch was used to extract the hard 
X-ray source position at R.A. (2000.0) = 17''25'"10' and Dec. 
- -36°17'18", with an uncertainty of 23" that includes a 
10" systematic error due to the instrument misalignment (see 
Fig. 1). This statistical error is of the same order as the system- 
atic error that comes from the image reconstruction, which is 
proportional to the source significance: ~ 36" for a S/N of 47 
(Gros et al. 2003). This position is compatible with the position 
obtained by Walter et al. (2004). 

Fig. 2 shows the XMM-Newton EPIC/MOSl image with 
the ISGRI error circle. A single X-ray source corresponds 
to the INTEGRAL position. Its position is R.A. (2000.0) = 
17''25'"11.4' and Dec. = -36°16'58.6" with an uncertainty of 
4". 

With the accurate source position obtained with EPIC, the 
Two Microns Afl-Sky Survey (2MASS) catalog (Cutri et al. 
2003) was searched for an infra-red counterpart. Only one 
2MASS source appears fully within the 4" EPIC error circle, 
2MASS J17251 139-3616575, located 1" away from the EPIC 
position. The source is not detected in the J-band with an upper 
limit of 14.2 magn at 95% confidence level, but appears in the 
H-band with 1 1.8 magn and the K^-band with 10.7 magn (see 
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Fig.2. IGR J17252-3616 EPIC/MOSl image with the 23" 
INTEGRALflSGRl error circle 
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Fig. 3. Infrared counterpart of IGR J 1 7252-36 1 6. Images taken 
from the 2MASS survey: J-band (left), H-band (middle) and K- 
band (right). The EPIC error circle of 4" is displayed. 2MASS 
J17251139-3616575 is only present in the H and K images. 



Fig. 3). The colours obtained when dereddening those infra- 
red magnitudes can either be interpreted as a close cool star or 
a distant hot star (Walter et al. 2005). 

4.2. Timing analysis 
4.2.1. Long-term variability 

The long term variability of IGR J17252-3616 was studied us- 
ing ISGRI data in two energy ranges: 20-60 keV and 60-150 
keV. Mosaic images were built for each three days revolution 
and light curves extracted for both energy bands. When the 
source was not detected, 3<t flux upper limits were calculated. 
Net exposures higher than 5 ks at the source location were con- 
sidered. Typical net exposures for one revolution vary from a 
few ks to several tens of ks. The source was never detected 
above 60 keV. The 3cr 60-150 keV average count rate upper 
limit is ~0.3 counts s"' (~7 mCrab) with the highest and lowest 
values being ~0.6 counts s ' (~14 mCrab) and ~0.6 counts s ' 
(~3.5 mCrab). The 20-60 keV light curve on the time scale of 
a spacecraft revolution is shown in Fig. 4. 

The source is detected in almost every single revolution. Its 
average 20-60 keV flux is 0.96+0.01 counts s"' (-6.4 mCrab), 
derived from a 20-60 keV mosaic image generated with data 



XMM 



52700 52800 52900 53000 53100 53200 53300 



Fig. 4. IGR J17252-3616 20-60 keV Hght curve with average 
count rate per revolution (~3 days). 3cr upper limits are indi- 
cated when the source was not detected. The XMM-Newton ob- 
servation time is shown. 



of the first three visibility periods for a total net exposure of 
1.7 Ms. The average flux of the object selecting only revolu- 
tions when the source is detected is 1.39+0.02 counts s"' (~9.3 
mCrab) as calculated from Fig. 4. When selecting revolutions 
for which only upper limits are available (rev. 46, 49, 50, 56, 
118 and 183) and building a mosaic image with a net exposure 
of 165 ks, the source is detected with a significance of 5cr, and 
the mean 20-60 keV flux reaches 0.19+0.04 counts s"' (-1.3 
mCrab); this is five times less than the average mean flux over 
all the revolutions and seven times less than the average flux 
over revolutions when the source was detected. Nevertheless, 
the source is persistent as, once cumulating enough exposure 
time, it remains detectable. 

On revolution time scales, the source flared up to 4 
counts s"' (-27 mCrab), or a factor of 10 stronger than the low- 
est detection of 0.4 counts s"' (-2.7 mCrab), on one occasion 
(see Fig. 4). IGR J17252-3616 varies by a factor -4 outside of 
this flare. 

Variability on time scales shorter than a single revolution 
was also investigated. Light curves based on flux per pointing 
were built considering the full INTEGRAL data set in the 20- 
60 keV energy range. The source is often not detected in single 
pointings. The average count rate per pointing when the source 
is detected with a significance higher than 4cr varies between 2 
and 4 counts s"' (-13-27 mCrab) outside the flares. Four flares 
are detected in total, three in the 2nd visibihty period and one 
in the 4th visibility period. The brightest flare occured at MJD 
52879 when the flux increased from -2 counts s"' to a peak of 
10 counts s"' (-67 mCrab) and decreased again in -1 day (see 
Fig. 5). The other flares reached a flux level of -6-8 counts s~' 
(-40-53 mCrab) and lasted less than one day. Considering the 
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Fig. 5. IGR J17252-3616 20-60 keV light curve with average 
count rates per pointing (~30 min). Around MJD 52879, we 
observe the brightest flare during the four visibility periods. Sc- 
upper limits are indicated when the source was not detected. 

observed count rates limit values of 0.2 and 10 counts s"', the 
source varies by a factor of more than 50 between different 
epochs. 

4.2.2. Pulse period 

Fig. 6 shows the XMM-Newton light curve. A periodic oscil- 
lation of ~ 400 s is clearly visible. A power spectral density 
distribution was created from the entire light curve in order to 
search for periodicity. A main peak was found around ~400 s 
with a harmonic at 200 s. 

Starting from the estimated period, the best period was 
searched for in the distribution when folding the light curve 
over a range of different periods. A significant peak was ob- 
tained, giving the best period of 414.8+0.5 s when fitted with 
a Gaussian (see Fig. 7 top). The uncertainty on the period 
was estimated using Eq. 14 of Horne & Baliunas (1986) and 
Eq. 2 of Hill et al. (2005). Lomb-Scargle periodograms were 
generated using Press & Rybicki (1989) fast method. The light 
curve was folded with the best period to obtain the pulse pro- 
file between 0.4-10 keV (see Fig.7 bottom). The XMM-Newton 
folded light curve is defined so that the minimum flux is set 
at phase corresponding to MJD 53085.55144. A first local 
maximum is observed at phase ~ 0.15 followed by a broad 
maximum at phase -0.6. The pulse fraction, defined as Pf - 
(/max - /min) /(/max + /min) wifli /^ax and /min being the maximum 
and minimum intensities of the folded light curve, respectively, 
reaches 68 ± 3%. 

An ISGRI 50 s binned light curve was made during the 
bright flare observed at MJD 52879 in two energy bands 20-40 
keV and 40-60 keV. Epoch folding the INTEGRAL data, start- 



15 - 




(Time 53085.55911 MJD) [sccl 



Fig. 6. IGR J17252-3616 EPIC/pn 10 s binned light curve. The 
gaps correspond to the periods with enhanced background ac- 
tivity that were discarded before starting the analysis. 

ing with the period found in the EPIC/pn light curve, returned 
a period of 413.7+0.3 s in the 20-40 keV band (see Fig. 8 top). 
The signal to noise of the 40-60 keV light curve is not signif- 
icant enough to detect the modulation. The pulse profile shape 
shows a broad peak without complex structures (see Fig. 8 bot- 
tom). The pulse fraction reaches 56+12% and is consistent with 
the value observed between 0.4-10 keV. 

The shape of the X- and y-ray pulse profiles show some dif- 
ferences. Due to the less significant detection in the ISGRI data, 
the number of phase bins is different in the pn and ISGRI folded 
light curves. The main peak in ISGRI seems broader than in pn. 
There is also a hint of a secondary peak in the pn data that is 
absent in the ISGRI data. The folded light curves cannot be 
compared in phase, because the phase of both pulse profiles 
do not correspond in time as the observations are separated by 
7 months. Notice that the 20-60 keV source intensity was three 
times brighter during the flare (MJD 52879) than during the 
XMM-Newton observation (~MJD 53085.6) (see Fig. 4). 

4.2.3. Orbital period 

Indications of an orbital period in the 20-60 keV light curve 
were searched for, based on flux per pointing. A Lomb-Scargle 
periodogram was generated between 4 and 28 days. The peri- 
ods with flaring activity were removed in the light curve before 
searching for any coherent modulation. A significant peak at 
9.72+0.09 days is visible and the folded light curve is shown in 
Fig. 9. A minimum flux consistent with an eclipse is visible at 
phase 0. It lasts between phases 0.97-1.1 that corresponds to a 
lapse of time of ~ 1 .26 days. The normalized intensity smoothly 
increases from phase 0.1 to 0.3 to reach a plateau that lasts 
~3.9 days. Afterwards, the intensity gradually decreases from 
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Fig. 7. top: distribution of trials periods for epoch fold- 
ing search on EPIC/pn data. The best value was fitted with 
a Gaussian. The best period of 414.8+0.5 s is indicated with 
a dashed-line. bottom: IGR J17252-3616 pulse phase folded 
light curve in the 0.4-10 keV energy band obtained with the 
best spin period. The zero epoch is MJD 53085.55144. 



Fig. 8. top: IGR J17252-3616 IBIS/ISGRI epoch folding 
search x~ distribution. The best period was fitted with a 
Gaussian and resulted in 413.7+0.3 s (see dashed-line). bot- 
tom: ISGRI folded light curve built with data from revolution 
106 (MJD 52877.4-52880.4). IGR J17252-3616 pulse phase 
folded light curve in the 20-40 keV energy band was obtained 
with the best spin period. The zero epoch is MJD 52877.50072. 



phase 0.7 to the start of the eclipse at phase 0.97. The XMM- 
Newton observation corresponds to phase ~0.7. 



4.3. Spectral analysis 

The X-ray spectral bins were grouped to have at least 100 
counts per channel (see Fig. 10). This allows the use of the;if^ 
statistic. The spectrum is strongly absorbed at energies lower and ISGRI spectra are used together 



than 4 keV. An emission line at 6.4 keV and an edge at 7 keV 
are clearly detected. The spectra were first fitted with a sim- 
ple absorbed power law and a blackbody model. In terms of 
reduced j^'^, these simple models fitted the XMM-Newton data 
well. In order to distinguish between the two models, the EPIC 
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Fig. 9. IGR J17252-3616 orbital folded light curve in the 20- 
60 keV energy band was obtained with the best period of 
9.72 d. The zero epoch is MJD 52671. The XMM-Newton ob- 
servation correspond to phase ~0.7. 

The ISGRI spectrum was extracted from revolution 106 
data when the source intensity was the strongest (see Fig. 5) 
and also from an average mosaic image of all visibility periods. 
The two spectra do not show any significant difference beside 
their normalisations. In the following spectral analysis, the rev- 
olution 106 source spectrum was used since it has the highest 
S/N. 

When fitting EPIC and ISGRI spectra together, an inter cal- 
ibration factor, CisGRi, was applied to take into account the dif- 
ferent source flux levels during the non simultaneous observa- 
tions. The simple absorbed power law and blackbody models 
failed to fit the data between 1-50 keV. The blackbody is too 
narrow and the power law requires a high-energy cutoff. 

A cutoff power law (CPL) was therefore chosen as a 
phenomenological model to fit the data. Moreover, a more 
physical model based on Comptonization (CTT) (Titarchuk 
1994) was also adopted as suggested by other observa- 
tions (Torrejon et al. 2004; Bodaghee et al. 2005; Maselti et al. 
2005; Walter et al. 2005). In Xspec, the CPL model is defined 
as cst"wabsCbbody+vphabs(cuto£fpl+ga)) and CTT as 
cst''wabs(bbody+vphabs(corapTT+ga)). The 6.4 keV fine 
was fitted with a Gaussian. An absorption model was used, 
where the iron abundance was left as a free parameter. The soft 
excess at low energies was modelled with a black body. The 
galactic absorption is also taken into account. All the spec- 
tral parameters errors were calculated at the 90% confidence 
level. Parameters resulting from the spectral fitting are listed in 
Table 2. 

The CPL model gives typical values of an accreting pul- 
sar with a flat powerlaw, F = 0.0 + 0.1, and an energy cut- 
off, £c = 8.2 + 0.4 keV. The CTT model also fits the data 




channel energy (keV) 



Fig. 10. Combined IGR J17252-3616 EPIC-nlSGRI averaged 
spectrum. All three EPIC spectra are displayed in the X-ray 
band: pn (top) and MOS[12] (bottom). The ISGRI spectra is in 
the hard X-ray band. The source is detected up to 50 keV. CTT 
model is used. 

well with characteristic values for the electron temperature of 
kT^ = 5.5 + 0.2 keV and an optical depth of t = 7.8 + 0.6. In 
both cases, the unabsorbed flux between 2-12 keV and 13-100 
keV is 2 and 5 10"'" ergs cm"^ s"', respectively. 

The line centroid at 6.40 + 0.01 keV is compatible with 
what is expected from cold iron. The line width is consistent 
with zero and was fixed for the spectral fits. The equivalent 
width of that line is 84 eV and was calculated considering the 
unabsorbed continuum. 

A soft excess is required below 3 keV and was represented 
by a blackbody absorbed by the galactic column density (with- 
out soft excess, ;\fVd.o.f.=415/376 and 420/376 for CPL and 
CTT models, respectively). In this case, the 1-2 keV soft X-ray 
flux is 1.6 or 2.0 10""^ ergscm^^ ^-i ((-pL or CTT). An alter- 
native model was also used to explain the origin of the soft 
X-ray excess based on the idea that the absorbing matter could 
only partially cover the X-ray emitting source (model pcfabs 
in Xspec). A covering factor of 0.995 +0.002 could also explain 
the soft excess with a good fit ofx^ld.o.f. = 406/378 for CTT 
model (41 1/378 for CPL). 

For phase-resolved spectroscopy with EPIC/pn, the pulse 
profile was cut into three phase intervals; 0.00-0.43 + 0.81- 
1.00 for the low intensity par;, 0.43-0.49 + 0.71-0.81 for the 
wings of the broad pulse; 0.49-0.71 for the pulse core. The ef- 
fective exposures for each phase resolved bin are 3383 s, 894 s 
and 1229 s, respectively. The phase-resolved spectra were ex- 
tracted and the channel bins were grouped to have at least 50 
counts per channel for the low intensity part and the pulse core, 
and 20 counts per channel for the pulse wings. All the result- 
ing spectra were fitted inside Xspec using the CTT model again 
combined with the ISGRI spectrum to better constrain the spec- 
tral fit. No significant variation of the spectral shape with the 
phase was found, excepting the change in the normalisation 
(see Table 3). 

Therefore, the three phase-resolved spectra were fitted with 
the Comptonization model and all the continuum parameters 
were fixed, except its normalisation, the column density, and 
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Table 2. Spectral analysis. The EPIC and ISGRI spectra were fitted together. Two models in Xspec were se- 
lected to fit the data. Model CPL is constant*wabs(bbody+vphabs(cutoffpl+gaussian)) and model CTT is 
constant*wabs(bbody+vphabs(compTT+gaussian)). The galactic absorption N^' was fixed. The errors are calculated at 
the 90% confidence level. 



Models 


Parameters 


Values 


Unit 


model CPL 


^ISGRl 


2.6 








1.5 


1022 gjjj-2 




soft 


0.5 


keV f fixed') 






1.6 


lO"''' ergs cm"^ s"' 






12 8+° '' 
1 5+"' 

-0.2 


10^2 cm-2 






Zo 




r 


0.02+!;" 

-0.10 








8.2+j;* 


keV 






6.400+^-!!?; 

" ""-0.013 


keV 






—yj.z. 


lO-^'phcm-^ 




unabs -FT-nkeV 


1.5 


10 '" ergs cm"^ s"' 




UTiabs mm \/ 
uiiciLiLi 1 IJ-lUUkeV 


5.1 


lO"' ergs cm^^ s ' 




r/d.o.f. 


401/376 




model CTT 


CiSGRl 


1.7 








1.5 


1022 cm"^ 




soft 


0.5 


keV ffixed'l 




^ l-2keV 


2.0 


1 n~ pro*; rrn~^ 




A^H 


15.3+!-' 

-1.0 


10-^ cm-^ 




Zpe 


1 4+0.3 
^•^-0.3 


Zo 




kTo 


0.1 


keV (fixed) 




kT, 


n c+0.6 

'•'^-0.5 


keV 




T 






^line 


5 401+0.005 


keV 






2 r" - 


10""' phcm"2 ^-1 




unabs F2_i2kev 


1.7 


10"'" ergs cm"2 s"' 




unabs Fi3_iookev 


4.8 


10"'" ergs cm"^ s"' 




^Vd.o.f. 


401/376 





Table 3. Phase-resolved spectral analysis. We used the CTT model. All continuum and normalisation parameters were left free. 
A few parameters were fixed as for the average spectrum (see Table 2); A^^'*' = 1.5 10^^ cm"^, kT^ofi - 0.5 keV, kTo = 0.1 keV. 
Only pn data were used for the phase-resolved spectral analysis. 



Parameters low intensity pulse wings pulse core Unit 



ClSGRl 


2.6 


1.4 


0.9 




/^l-2kcV 


1.4 


4.2 


2.8 


10 ''' ergs cm - s ' 


A'h 


15 2+2-" 


13 7+" 


15 r' ** 


10^2 cm"2 


Zpe 


1 9+0.6 
^•^-0.5 


2 0+'-^ 


1 6+"-' 

1-O_0.4 


Zo 


kT, 


5 5+"-2 




5 5+"-2 


keV 


T 


8.0-g 


8.7!|i 


7 S+0.8 
'•°-0.7 




^line 


6.40!jJ:«i 


6 42+"-"" 

"•^^-0.03 


6 40+" "'' 
"•^"-0.05 


keV 


Fv 
^ line 


9 9+0.3 
^■^-0.3 




1 4+"-^ 

^- -0.6 


10 -* phcm-2 s ' 


unabs F2-i2kcV 


1.2 


2.0 


3.4 


10"'" ergscm"2 s"' 


X^I±o.i. 


145/160 


208/217 


180/197 





the soft excess and the line normaUsations. The results are sum- 
marized in Table 4. No significant variation of the line flux nor 
of the absorbing column density was observed. The unabsorbed 
2-10 keV flux changes according to the pulse phase. The varia- 
tions of the 1-2 keV flux are not significant once the errors are 
considered. 



5. Discussion 

5.1. X-ray counterpart 

When IGR J17252-3616 was discovered in the y-rays, it was 
first labelled as a new source since no other known sources 
lay within the ISGRI error circle of 2'. The nearest source 
was EXO 1722-363 at a distance of 12'. The EXOSAT source 
was considered to be too far to be the X-ray counterpart of 
IGR J17252-3616. With the improved ISGRI and EPIC po- 
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Table 4. Phase-resolved spectroscopy with all the spectral modeling parameters of model CTT except normalisations equal to 
the values of the average spectrum (see Table 2). Only pn data were used for the phase-resolved spectral analysis. 



phase 




10-"* phcm-2 s-i 


Fl-2 kcV 

10"'* ergs cm"- s"' 


unabsorbed Fo-io kcV 
10"'" ergscm"^ s"' 


x'/d.o.f 


low intensity 
pulse wings 
pulse core 


15 3+'" 
15.8+i;^ 
15 5+°-' 

1 J.J Q J 


9 1 +0.4 
^■^-0.3 
1 5+"' 

1 6+""^ 

^■"-0.7 


^■•^-1.2 

6.8!°-^ 
2-8!?i 


0.9 
1.4 
2.4 


147/164 
214/221 
181/201 




2 5 10 20 

channel energy (keV) 



Fig. 11. EPIC/pn phase-resolved spectra. Only the pn spectra 
of the three phase intervals are displayed in the X-ray band: 
low intensity (bottom), pulse wings (middle) and pulse core 
(top). The ISGRI spectrum in the hard X-ray band is not phase 
resolved. 

sition, the source still lies 12' away from the given EXOSAT 
position with an accuracy of 9' (Warwick et al. 1988). 

Tawara et al. (1989) also gave a new position within a box 
of 2° X 4° that is 20' away from the EPIC position. They dis- 
covered a pulsation of 413.9 s. The pulse profile shows a single 
peak that is independent of the energy and a high pulse frac- 
tion of 80%. The spectrum is heavily absorbed at low energy 
and shows a very hard shape plus an emission line at 6.2+0.5 
keV. Even if the position given by all these studies is quite far 
from the accurate IGR J 17252-36 16 position, the informations 
given by Tawara et al. (1989) strongly confirm that the same 
source is seen since the same features are seen in this study. 

From archival data of RXTE/PCA taken between 1999 
and 2003, Markwardt & Swank (2003) reported an orbital pe- 
riod of 9.737+0.004 days, consistent with the lower limit 
given by Takeuchi et al. (1990). This result with RXTE was 
confirmed by Corbet et al. (2005). Their best orbital pe- 
riod is 9.741+0.004 d. They also find a pulse period of 
413.88553+0.00001 s, equivalent to the previous studies and 
with the measurement presented here. 

5.2. Temporal modulations 

In this paper, periodic or episodic time variations were searched 
for on time scales of seconds to days. First, the two years long 
ISGRI 20-60 keV light curve indicates that the source is per- 
sistent. However, variations in the y-ray count rates larger than 



50 were observed for this system. Four flares were detected on 
a time scale of hours. 

An orbital period of 9.72+0.09 days was derived that is 
consistent with the one derived by Corbet el al. (2005). The 
eclipse duration is a bit lower in the ISGRI data than in the 
RXTE data (~1.3 days and ~1.7 days respectively). Both or- 
bital profile at 2-10 keV and 20-60 keV are consistent. The 
progressive decrease of the intensity before entering the eclipse 
between phase 0.7 and 0.97 can be explained by the changing 
hydrogenic column density that increases from ~ lO^"* cm"^ to 
~ 10""* cm"^. This larger amount of matter along the line of 
sight implies that a larger part of the flux is scattered. 

Similar pulsations were found in EPIC/pn and ISGRI data: 
414.8+0.5 s and 413.7+0.3 s, respectively. The EPIC and 
ISGRI pulsations are consistent with the one detected by 
Tawara etal. (1989) and Corbet et al. (2005). Therefore, the 
spin period did not change significantly in the last 17 years. 
Takeuchi et al. (1990) also searched for pulse period changes 
and did not find any significant spin variation in Ginga data. 
The neutron star seems to not accrete kinetic momentum, 
which suggests accretion from a stellar wind. 

With such spin pulsation and orbital period, the source 
is situated in the underfilled Roche-Lobe supergiant region 
in the P^pm vs forb Corbet diagram (Corbet 1986). Together 
with the fact that the source is persistent and shows non- 
periodic flares with different intensities, this strongly suggests 
that IGR J 17252-36 16 is a high mass system fed by stellar 
wind with the primary star being a supergiant. 

5.3. Spectra 

The EPIC and ISGRI spectra that could be well fitted with a 
flat power law (F ~ 0) and an energy cutoff at Ec = 8.2 keV 
are typical of X-ray pulsars (White et al. 1995). Three spectral 
features can give more clues about the physics of this object: 
the apparent soft excess, the huge hydrogen column density and 
the cold Fe Ka line. 

A soft excess is detected below 2 keV. Such soft X-ray ex- 
cess has been observed in other HMXB (Hickox et al. 2004). 
It could originate in X-ray scattering or partial ionization in 
the stellar wind (White el al. 1995). The spectrum at soft X- 
rays could also be explained by a partial covering of the X-ray 
source by the absorbing matter. In the latter case, pulsations at 
soft X-rays would be expected. The 0.4-2 keV light curve was 
extracted and folded using the same ephemeris and period used 
to fold the 0.4-10 keV light curve. Then, it was compared to 
a constant model withj^^^/d.o.f. ~ 57/50 and the 0.4-10 keV 
pulse profile with;^'-/d.o.f . ~ 54/50. Both models are compati- 
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ble with the 0.4-2 keV folded light curve. However, the number 
of events gathered below 2 keV is too poor to decide if the soft 
excess is pulsating or not. 

The absorbing column density was estimated at 13-15 10^^ 
atoms cm"^ (model CPL-CTT) that is ten times larger than ex- 
pected on the line of sight A^h = 1-5 10^^ atoms cm"^. This 
high absorption could be explained by the stellar wind ex- 
pelled by the primary star that surrounds the neutron star. 
During the XMM-Newton observation, the column density did 
not vary with the pulse. The fact that the amount of matter 
did not change noticeably during the 10 ks XMM-Newton ob- 
servation could indicate that the surrounding matter is stable 
on this time scale. However, other values of the column den- 
sity were reported at different epochs by Tawara et al. (1989), 
Takeuchi et al. (1990) and Corbet etal. (2005). The column 
density has been observed to increase up to lO^'* atoms cm"^ in 
the past. This could be related to the neutron star moving along 
the orbit (Tawara etal. 1989; Corbet et al. 2005). Therefore, 
the matter does not homogeneously surround the binary sys- 
tem. Instabilities in the stellar wind of the companion star could 
also be responsible for this evolving absorption. However, in 
Corbet et al. (2005), the column density seems to be linked to 
the orbital phase where the highest values correspond to the exit 
of the eclipse and the lowest ones when the source is in front of 
the companion star. The absorption is therefore intrinsic to the 
binary system. 

The 6.4 keV Fe Ko' line is detected. The Fe Ka line flux 
does not vary with the pulse phase (see Table 4). Therefore, 
either the matter responsible for the fluorescence is homoge- 
neously distributed around the source or the thickness of the 
shell emitting the fluorescence is larger than 1 .2 10^ km around 
the accreting system. The Fe Ka line energy is 6.401;^Q Qgg keV 
(CTT model) and corresponds to iron that is at most 12 times 
ionized (House 1969). An upper limit of the ionization param- 
eter can be estimated as log(L/nR^) < 1, where L is the lumi- 
nosity, n the gas density, and R is the distance from the ioniz- 
ing continuum source to the inner shell surface (Kallman et al. 
2004). Since nR ~ A^h, the distance of the fluorescence source 
from the X-ray source is larger than 10^ km. The values de- 
rived for the hydrogen column density, the equivalent width, 
and the estimated over abundance of iron are compatible with a 
spherical distribution of matter around the source (Matt 2002). 

Considering a typical luminosity of an active accretion- 
powered pulsars of 1.2 10^^ ergs s"' (Bildsten et al. 1997), the 
distance of IGR J17252-3616 can be estimated as 7 kpc, close 
to the galactic centre. 

6. Conclusions 

IGR J 17252-36 16 is the hard X-ray counterpart of 
EXO 1722-363, and the most accurate source position 
to date has been provided here. IGR J17252-3616 has been 
monitored by INTEGRAL during two years for a total exposure 
of 6.5 Ms and XMM-Newton performed a follow-up obser- 
vation of 3 hours. The source is persistent with an average 
20-60 keV flux of 6.4 mCrab. Four flares lasting less than 
one day were detected by INTEGRAL. The source's count 
rate varies by a factor larger than 50 on such timescales. A 



pulsation has been detected in both EPIC and IBIS/ISGRI 
data of 414.8+0.5 s and 413.7+0.3 s, respectively. There is 
no evidence of spin period variation. An orbital period of 
9.72+0.09 d is also found in IBIS/ISGRI data. The spectral 
shape is typical for an accreting pulsar except that a huge 
intrinsic absorption and a cold iron fluorescence line are 
detected. The absorbing column density and cold iron line 
do not vary with the pulse period. The absorbing fluorescent 
material is distributed around the neutron star in a shell com- 
parable in size with the orbital radius. With the accurate X-ray 
position, we provide a likely infra-red counterpart within the 
X-ray error box: 2MASS J17251 139-3616575. The observed 
features of the source suggest that it is a wind-fed accreting 
pulsar. This object is a new member of the growing family 
of heavily-obscured HMXB systems that have been recently 
discovered with INTEGRAL. The source is located ~7 kpc 
away, near the Galactic Centre. 
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